Inhomogeneous freeze-out in relativistic heavy-ion collisions 
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A QCD phase transition may reflect in a inhomogeneous decoupling surface of hadrons produced 
in relativistic heavy-ion collisions. We show that due to the non-linear dependence of the particle 
densities on the temperature and baryon-chemical potential such inhomogeneities should be visible 
even in the integrated, inclusive abundances. We analyze experimental data from Pb-(-Pb collisions 
at CERN-SPS and Au-I-Au collisions at BNL-RHIC to determine the amplitude of inhomogeneities. 
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INTRODUCTION 

Collisions of heavy nuclei at relativistic energies pro- 
duce very hot and baryon-dense QCD matter [1]. The 
entropy per net baryon in the central region increases 
monotonically with beam energy, i.e. the maximum tem- 
perature increases, while the baryon-chemical potential 
decreases. Hence, by varying the energy as well as the 
mass number of the nuclei etc., one can explore various 
regimes of the phase diagram of QCD. In particular, it 
is expected that at sufficiently high energies a transient 
state of deconfined matter with broken Z{3) center sym- 
metry and/or with (approximately) restored chiral sym- 
metry is produced. 

Lattice QCD simulations [2] indicate that a second- 
order critical point exists, which was predicted by effec- 
tive chiral Lagrangians [3]; present estimates locate it 
at T « 160 MeV, fin ~ 360 MeV. This point, where 
the (T-field is massless, is commonly assumed to be the 
endpoint of a line of first-order phase transitions in the 
{fJ-B,T) plane. 

There is an ongoing experimental effort to detect the 
line of first-order phase transitions at high baryon den- 
sity, and its possible endpoint, in heavy-ion collisions. 
It is hoped that by varying the beam energy, for ex- 
ample, one can "switch" between the regimes of first- 
order phase transition and cross over, respectively. If the 
particles decouple shortly after the expansion trajectory 
crosses the line of first order transitions one may expect 
a rather inhomogeneous (energy-) density distribution on 
the freeze-out surface [4, 5] (similar, say, to the CMB pho- 
ton decoupling surface observed by WMAP [6]). On the 
other hand, if the low-temperature and high-temperature 
regimes are smoothly connected, pressure gradients tend 
to wash out density inhomogeneities. Similarly, in the ab- 
sence of phase-transition induced non-equilibrium effects, 
the predicted initial-state density inhomogeneities [7, 8] 
should be strongly damped. 

Here, we investigate the properties of an inhomoge- 
neous fireball at (chemical) decoupling. Note that if the 
scale of these inhomogeneities is much smaller than the 
decoupling volume then they can not be resolved indi- 
vidually, nor will they give rise to large event-by-event 



fiuctuations. Because of the nonlinear dependence of the 
hadron densities on T and /xs, they should nevertheless 
refiect in the average abundances. Our goal is to check 
whether the experimental data shows any signs of inho- 
mogeneities on the freeze-out surface. 

Our basic assumption is that as the fireball expands 
and cools, at some stage the abundances of hadrons 
"freeze" , keeping memory of the last instant of chemical 
equilibrium. This stage is refered to as chemical freeze- 
out. By definition, only processes that conserve particle 
number for each species individually, or decays of unsta- 
ble particles may occur later on. 

The simplest model is to treat the gas of hadrons 
within the grand canonical ensemble, assuming a homo- 
geneous decoupling volume. The abundances are then 
determined by two parameters, the temperature T and 
the baryonic chemical potential fis', the chemical poten- 
tial for strangeness is fixed by the condition for overall 
strangeness neutrality, and other conserved charges shall 
be neglected. Fits of hadronic ratios were performed ex- 
tensively [9, 10] within this model, sometimes also includ- 
ing a strangeness (7s) or light quark (7^) supression fac- 
tor [11, 12] or interactions with the chiral condensate [13]. 

The purpose of the present paper is to analyze the ex- 
perimental data on relative abundances of hadrons with 
respect to the presence of inhomogeneities on the decou- 
pling surface. To that end, we propose a very simple 
and rather schematic extension of the common grand 
canonical freeze-out model, i.e. a superposition of such en- 
sembles with different temperatures and baryon-chemical 
potentials. Each ensemble is supposed to describe the 
local freeze-out on the scale of the correlation length 
~ 1/T '--^ 1 — 2 fm. Even if freeze-out occurs near the 
critical point, the correlation length of the chiral conden- 
sate is bound from above by finite size and finite time 
effects, effectively resulting in similar numbers [14]. On 
the other hand, for small chemical potential, far from 
the region where the a-field is critical, the relevant scale 
might be set by the correlation length for Polyakov loops, 
which is of comparable magnitude [15]. Classical nucle- 
ation theory for strong first-order phase transitions pre- 
dicts even larger "bubbles" [16] but is unlikely to apply to 
small, rapidly expanding systems encountered in heavy- 
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ion collisions [5, 17]. Another (classical) model for the 
formation of small droplets in rapidly expanding QCD 
matter has been introduced in [18]. 

The entire decoupling surface contains many such "do- 
mains" , even if a cut on mid-rapidity is performed. We 
therefore expect that the distributions of temperature 
and chemical potential are approximately Gaussian. Be- 
sides simplicity, another goal of the present analysis is to 
avoid reference to a particular dynamical model for the 
formation or for the distribution of density pcrtmbations. 
In fact, we presently aim merely at checking whether any 
statistically significant signal for the presence of inhomo- 
gencities is found in the data. If so, more sophisticated 
dynamical models could be employed in the future to 
understand the evolution of inhomogeneities from their 
possible formation in a phase transition until decoupling. 

Rate equations for nuclear fusion and dissociation pro- 
cesses (and neutron diffusion) have been used for inho- 
mogencous big bang nucleosynthesis in the early uni- 
verse [19]. Similarly, hadronic cascade models could be 
used for heavy- ion reactions [20] . This would remove ref- 
erence to the grand-canonical ensemble and to a thin de- 
coupling surface in space-time. In fact, hadronic binary 
rescattering models do predict a rather thick freeze-out 
layer [20, 21], where matter expands non-ideally. On the 
other hand, the steep drop of multi-particle collision rates 
with temperature should narrow the freeze-out again [22] . 
In either case, we do not expect a strong energy depen- 
dence of the width of freeze-out (see also [23]). 

At chemical freeze-out, matter is in a state of expan- 
sion. However, such flow effects do not affect the rela- 
tive abundances of the particles (in full phase space) if 
their densities are homogeneous throughout the decou- 
pling volume. The total number of particles of species i, 
integrated over a solid angle of 47r, is given by an integral 
of the current TVf = PiU^, with the four- velocity of 
the expanding fluid, over a given freeze-out hypersurface 



and so can be written in the form 



u^da,. 



(1) 



The second factor on the r.h.s. is nothing but the three- 
volume V3 of the decoupling hypersurface as seen by the 
observer. This volume is common to all species and drops 
out of multiplicity ratios: Ni/Nj = p^°/p^°. It is clear 
that the argument holds even when cuts in momentum 
space are performed, provided that the differential distri- 
butions of all particles do not depend on that particular 
momentum-space variable (for example, rapidity cuts for 
boostinvariant expansion [24]). 

When the intensive variables T and pB vary, then the 
integration measure u- du)/Vz will, in general, depend 
on the assumed distribution and amplitude of inhomo- 
geneities, as well as on the hydrodynamic flow profile etc. 
Nevertheless, it is still the same for all particle species 



^Ju-da — > J dTdpB P{T, pb) 



(2) 



with P{T, pb) some distribution for T and Pb- For sim- 
plicity, and for lack of an obvious motivation for assum- 
ing otherwise, we shall take P{T,pb) to factorize into 
a distribution for T, times one for pB- These distribu- 
tions could, in principle, be obtained from the real-time 
evolution of the phase transition [4, 5] . 



THE MODEL 

With these qualifications aside, we now proceed to in- 
troduce our model and to analyze the available data from 
heavy-ion collisions at high energies. The hadron abun- 
dances are determined by four parameters: the arith- 
metic means of the temperatures and chemical poten- 
tials of all domains, T and Pb, and the widths of their 
Gaussian distributions, ST and SpB- Then, the average 
density of species i is computed as 



Pi iT,PB,ST,5pB) = 



(3) 



J dT P{T;T,6T) j dpB P{pb; Pb^^I^b) Pi{T, pb) , 
-00 

with Pi(r, Pb) the actual "local" density of species i, and 
with P{x;x,5x) ^ exp[— (a; — S)'^/2 iJo;^] the distribu- 
tion of temperatures and chemical potentials within the 
decoupling three-volume (the proportionality constants 
normalize the distributions over the intervals where they 
arc defined). In the limit 5T. Sps ^ the Gaussian 
distributions are replaced by J-functions and the con- 
ventional homogeneous freeze-out scenario is recovered: 
Pj (T, /Is, 0,0) = /9,;(T./Ig). In other words, in that limit 
the average densities are uniquely determined by the first 
moments of T and Pb- 

For the present analysis we compute the densities 
pi{T,pB) in the ideal gas approximation, supplemented 
by an "excluded volume" correction: 



Pi{T,pB) = - 



id— gas 



{T, pb) 



+v^I: 



id- 



-gas 



(4) 



This schematic correction models repulsive interactions 
among the hadrons at high densities. Vi denotes the 
volume occupied by a hadron of species i; we employ 
V = |7ri?o^ with Rq = 0.3 fm for all species [25]. There- 
fore, for the homogeneous model the denominator in (4) 
drops out of multiplicity ratios. This is not the case for 
an inhomogeneous decoupling surface, where the distri- 
butions of various species diflPer. 

The densities of strange particles depend also on the 
strangeness-chemical potential ps, which we determine 
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by imposing local strangeness neutrality. Strictly speak- 
ing, this condition is not mandatory, of course. The ef- 
fect of independent fluctuations of /xg should be looked 
at in the future, in particular for collisions at low and 
intermediate energies {\/snn^ 15 GeV). This may help 
to reproduce the A to p ratio, which was found to be 
larger than one [26] and the X+/7r+ enhancement around 
El^^/A = 30 GeV [27]. For all fits over the full solid an- 
gle, we fixed the isospin chemical potential by equating 
the total charge in the initial and final states; for the 
mid-rapidity fits at high energies, we fixed /x/ = 0. 

To illustrate the eSect of inhomogeneities on the distri- 
butions of various hadrons within the decoupling volume 
we introduce 

Di{T;T,JlB,ST,S^iB) = P{T;T,ST) 

OO 

— , (5) 



Pi{T,JlB,ST,dfiB) 
Di{fiB;T,JlB,ST,6fiB) = P{ijlb;1j.b^^1^b) 

oo 

fdTP{T;T,5T) Pi{T,i^B) 

w 



Pt iT,PB,^T,SnB) 



(6) 



Di{T), for example, is the probability that a particle of 
type i was emitted from a domain of temperature T. 
The main contribution to the integrals in (3) is not from 
T and JIb since hot spots shine brighter than "voids". 
Rather, they are dominated by the stationary points of 
the distributions defined in eqs. (5,6) above. Hence, the 
average emission temperature {T)i and baryon-chemical 
potential {fiB)i in general depend on the particle species 
i, unless ST = Sp,B = 0. They can be evaluated as 



{T)i 



{l^B)i 



CSC 

-I 



dTT Di{T;T,TlB,ST,S^iB) 



dfiB Mb Di{iiB; T,Pb^ ^T, 6iib) ■ (7) 



Hence, masslcss particles arc not very sensitive to small 
fiuctuations in temperature. They are typically emitted 
from regions with temperature k, T, up to quadratic cor- 
rections in ST/T. In other words, the "particle emission 
distribution" D{T) is shifted (and skewed) only slightly 
from the "temperature distribution" P{T). This, of 
course, is simply due to the p ^ power-law form of 
the local density of massless particles. The convolution 
of a (narrow) Gaussian with a power-law does not lead 
to a large shift of the peak. 

Next, we consider the more interesting case of mas- 
sive particles with chemical potential equal to p.B (for 
anti-baryons replace JI^ — > etc.). We assume that 

exp[(/XB — m)/T] <C 1 and m/T ^ 1 such that quantum 
statistical and relativistic corrections can be neglected, 
but allow that x = {m/T){5T'^ /T^) =0(1). The integral 
over /xb is then straightforward, leading to the particle 
distribution 



D{T) ^ 

exp[^log(T/T) + ^ 



2T 



[t-tY 

2 6T'^ 



Again, we look for the stationary point of the exponen- 
tial. To leading order in the Gaussian width. 



= 1 + 



5T / m 
ST^ ( m 

T 



Mb 
T 



0(1) ~ l + x 



(11) 



In the second step, we have used the fact that in the 
Boltzmann limit Hb/T is of order 1 (i.e., not parametri- 
cally large). Hence, for massive particles the distribution 
is shifted by a large amount x = 0(1). Their emission is 
dominated by the tails of the Gaussian distribution or, 
in physical terms, by rare "hot spots" . 

To estimate the increase of the average density relative 
to the homogeneous case we plug the expression for T* 
into D{T) from eq. (10), which gives 



Physically, this means that for non-zero widths of the 
temperature and chemical potential distributions the 
freeze-out volume is not perfectly "stirred", in that the 
relative concentrations of the particles vary. 

For some limiting cases one can estimate the effect ana- 
lytically. Consider first massless particles without chem- 
ical potential (w direct pions): 



Z?(T;r,<5r)-expl 31og = 



T {T-T)' 



2ST^ 



(8) 



To leading order in 6T/T the stationary point of the ex- 
ponential is 

T 



= l + S +.. 



(9) 



■p{T,5TjiB,5iJLB) 
p{T,Jib) 

f exp 



1 SfXB"^ 

(1 + xY 2T' 



^ ~ Mb \(i2) 
1 + x T ' 



The non-exponential prefactor / arises due to the chang- 
ing width of the integration measure when the temper- 
ature is inhomogeneous and can not be estimated by a 
saddle-point integration (/ = 1 for dT = 0). Regardless, 
the main issue here is that this ratio increases exponen- 
tially with (SfiB/T)'^. Hence, we conclude that small fluc- 
tuations in the chemical potential are sufficient to raise 
the number of heavy particles significantly. Note that 
what matters is the magnitude of Sp.B relative to T, not 
Mb- 
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Similarly, temperature inhomogeneities also increase 
the density exponentially. That is, for 6T 0, x/{l + 
x) « X. However, the growth saturates when a; 1. The 
enhancement factor from a; = to a; ^ 1 is very large, 
~ Vexp((m-7lB)/T). 

DATA ANALYSIS 

To determine the four parameters of the model we min- 
imize 

^2^J2{rr-rr'''')'M (13) 

i 

in the space of T, Jl^, 6T, and 6^b- That is, we obtain 

least-square estimates for the parameters, assuming that 
they are independent. In (13), r"^^^ and r^°d.^^ denote 
the experimentally measured and the calculated particle 
ratios, respectively, and cr| is set by the uncertainty of the 
measurement. Wherever available, we sum systematic 
and statistical errors in quadrature. 

The data used in our analysis are the particle multi- 
plicities measured by the NA49 collaboration for central 
Pb-hPb collisions at beam energy i^LabM = 20, 30, 40, 80 
and 158 GeV [27], and those measured by STAR for cen- 
tral Au+Au collisions at BNL-RHIC, ref. [28] {^/snn = 
130 GeV, compiled in [29]) and ref. [30] (200 GeV). At 
RHIC energies, we analyze; the miclrapidity data; at top 
SPS energy, both, midrapidity and An data. At all other 
energies, we restrict ourselves to the An solid angle data 
by NA49 in order to avoid biases arising from differing 
acceptance windows of various experiments. Further- 
more, our checks showed that the fit results can depend 
somewhat on the actual selection of experimental ratios. 
Hence, where possible, we have opted for the least bias 
by choosing rf ^ = N-^'^/N^^'p, i.e. the multiphcity of 
species i relative to that of pions. This represents the 
maximal set of independent data points, as it is equiv- 
alent to fitting absolute multiplicities with an additional 
overall three-volume parameter, Ni = V^pi- 

Specifically, at E^Lab/^ = 20, 30, and 80 GeV the mul- 
tipHcities of 7r+, 7r~, K~ , B — B, A, A, and (f) are 
available. For the (in-)homogeneous model, this leaves 
five (three) degrees of freedom. At 40 GeV, we can add 
the and fi-|-f2. The data sets for top SPS energies in- 
clude yet a few more species: p, p (only midrapidity). Kg 

(only 471"), and fi, seperately. For RHIC-130, we 
fitted to the K+/K-, p/p, A/A, fi/il, R-fir-, 

Kg/n^, p/n^, A/tt~ , Kq/tt~ , (ji/n^ . S^/tt^ and Vl/tt~ 
ratios. Finally, at RHIC-200 the K+ / R- , p/p, Tl/n, 
K~ /■K~, p/tt", A/tt", A/tt", S~/7r~, !E!+/7r~, O/tt", 
4'/K~ and Kq/K~ ratios were used. The first three ra- 
tios are close to unity and essentially just set the chemical 
potentials to zero; they do not help to fix T, 5T and SfjiB- 
Where appropriate, feeding from strong and electro- 
magnetic decays has been included in rf^"'^^'- by replacing 



Pi Pi + Bij pj . The implicit sum over j ^ i runs over 
all unstable hadron species, with Bij the branching ratio 
for the decay j ^ i, which were taken from [31]. Prom 
all the resonances listed by the Particle Data Group [31], 
mesons up to a mass of 1.5 GeV and baryons up to a mass 
of 2 GeV were included, respectively. The finite widths 
of the resonances were not taken into account, and un- 
known branching ratios were excluded from the feeding. 
These details are irrelevant for the qualitative behavior 
of ST and SfiB but do, of course, matter for quantita- 
tive results. For SPS energies, feeding from weak decays 
was included only for the S — > A and O — > A channels 
(and the corresponding decays of the anti-baryons) at 
-S'LabM = 20 and 30 GeV. All other experimental yields 
were already corrected by NA49 for feed- down from weak 
decays [27]. At RHIC energies, feeding from weak de- 
cays has to be taken into account for a variety of particle 
species, c.f. [28, 30]. 

Technically, our fits use a four-dimensional lookup ta- 
ble for Pi with 1 MeV steps in T and 6T and 10 MeV 
steps in /x and 6p for the inhomogeneous fits, and 1 MeV 
steps for both T and /i for the homogeneous fits. The 
finite grid of course limits our ability to determe the best 
fit. However, the effect was checked to be small, so that 
the accuracy should be sufRcient to investigate the qual- 
itative behavior of inhomogeneities, and whether they 
can improve the agreement with the experimental data 
significantly. 



RESULTS 




FIG. 1: X^/dof versus ^sjvjv for the homogeneous {5T = 
Sp = 0, squares) and the inhomogeneous fit {5T and Sp, free 
parameters, circles). The full and dashed lines are meant to 
guide the eye. Furthermore, the /dof corresponding to the 
95.4% confidence interval is shown by the dotted line. 

Fig. 1 shows the minimal per degree of freedom 
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FIG. 2: X^ /dof contours in the ST, SjjiB plane for top SPS 
energy, BLab = 158 GeV. The other two parameters (T, Jib) 
are allowed to vary freely. The jdof minimum is indicated 
by the cross. 



(taken as the number of data points minus the num- 
ber of parameters) for the homogeneous and the inho- 
mogeneous approach, respectively. At £^Lab/^ = 20 GeV 
and at RHIC energies, ^ jdof is similar for both mod- 
els. Thus, the inhomogencous model docs not provide 
a statistically significant improvement of the description 
of the measured particle ratios. Hence, the assumption 
of a nearly homogeneous decoupling surface can not be 
rejected for low SPS and RHIC energies. 

On the other hand, for intermediate SPS energies, 
-E'Lab/^ — 30 — 160 GeV, x^/rfo/ is considerably smaller 
for the inhomogencous freeze-out surface than for the 
homogeneous case, which is far outside the 95.4% con- 
fidence interval [32]. This indicates that the parame- 
ters of the inhomogencous model are well determined 
and have reasonably small error bars. It is worth not- 
ing that, in general, the improvement is not driven by 
one single species; rather, the inhomogencous model de- 
scribes nearly all multiplicities better than a homoge- 
neous decoupling surface [33]. We also stress that our 
fits with (5T = = reproduce results from the liter- 
ature [10, 11, 30] if the same input ratios are selected. 

To illustrate the significance of inhomogcneitics differ- 
ently, we show contours of x^/do/ in the plane of 
in figs. 2 and 3. Here, T and /xb were allowed to vary 
freely such as to minimize at each point. Fig. 2 shows 
that is relatively flat along the b\XB direction, while 
is determined more accurately and is clearly non-zero. 
In general we find that there is little correlation between 
bT and b\XB and that about the minimum, ^ is rather 
flat in b\iB direction. 

On the other hand, Fig. 3 shows that at RHIC energy, 
X^ is very flat in both directions. With the present data 
points, a homogeneous freeze-out model appears to be 
a reasonable approximation at high energies. We have 
traced the origin of the very flat ^{5T) at RHIC en- 
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FIG. 3: Same as fig. 2 for RHIC energy {^snn = 200 GeV). 

ergy to a strong anti-correlation between 6T and T. This 

degeneracy appears, in the regime of small chemical po- 
tentials, when the data is dominated by particles of simi- 
lar mass (including contributions from resonance decays). 
There is then no distinguishing feature which would fix 
both ST and T independently, i.e. larger ST can be traded 
for smaller T and vice versa. Eq. (12) indicates that the 
model requires several species whose densities are domi- 
nated by significantly different mass scales. 

We now proceed to analyze the energy dependence of 
the freeze-out conditions. As already mentioned above, 
the fit parameters T, JIq, ST and SjJs do not acquire a 
direct physical meaning; for example, T is simply the 
arithmetic mean of the temperature within the entire 
volume, but particle production is dominated by "hot 
spots". Similarly, the RMS variation of the tempera- 
ture within the decoupling volume, ST, is not a direct 
measure for the different emission temperatures of vari- 
ous particles, since the latter are given by a convolution 
of the Gaussian temperature distribution with the re- 
spective particle densities pi, and so depend also on the 
particle masses, the chemical potentials etc. 

Hence, rather than discussing the energy dependence 
of the above technical parameters, we instead focus on 
the particle emission temperatures (T)j and baryon- 
chcmical potentials introduced in (7). For each en- 

ergy, we define an average particle emission temperature 
by summing over all species i {after resonance decays): 



(r) = 



(14) 



and similarly for (/Us). 

Along the same lines, we can determine the variance 
of the emission temperatures (and chemical potentials) 
of various particle species as 



S{T) 



{{Th-{T)y 



(15) 
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FIG. 4; Average particle emission temperature (T) and chem- 
ical potential (/is) at -ELabM = 20, 30, 40, 80, 158 GeV 
(CERN-SPS) and ^sjvjv = 130 and 200 GeV (BNL-RHIC), 
from right to left (symbols). The "error bars" indicate the 
RMS deviations over particle species (i.e. 5{T} and 5{fiB))- 
For SPS-20 and RHIC, these are depicted by dotted lines 
since the inhomogeneities are not statistically significant. 



We repeat that all of these quantities are entirely deter- 
mined by the parameters T, Jlq, ST and SfXs charac- 
terizing the Gaussian distributions of temperature and 
baryon-chemical potential; they do not represent addi- 
tional fit parameters. In particular, when ST = Sfis — 0, 
the homogeneous model is recovered: (T) = T^", S{T) = 
etc. 

Fig. 4 depicts the energy dependence of the average 
freeze-out temperature and chemical potential, as defined 
in eq. (14). The general trend is the same as in the ho- 
mogeneous model: (T) increases with energy, while {(J-b) 
decreases. Their values agree to within 15% with those 
for a homogeneous system obtained previously by oth- 
ers (cf. for example fig. 3 in Redlich et al. or fig. 27 in 
Braun-Munzinger et aZ. [10]). The inhomogeneous model, 
however, predicts sizeable variations of the emission tem- 
peratures of different particle species: this is indicated in 
fig. 4 by the "error bars", which correspond to S{T) and 
5{^b) as defined in eq. (15). They were obtained from 
the best fit at each energy, i.e. using the parameters cor- 
responding to the lowest x^/o?o/ from the inhomogeneous 
model. We repeat, however, that at the lowest and high- 
est energies the inhomogeneous fit has no statistical sig- 
nificance, in that /dof is only marginally better than 
for the homogeneous model. Hence, at these energies, the 
resulting S{T) and S{iib) (dotted lines) should be taken 
with ~ 100% error. 

At intermediate and high SPS energies, however, the 
variance of the temperature can be determined reason- 
ably well and is significantly larger than zero. S{T) ap- 
pears to be rather large already at i^LabM ^ 30 GeV, 
and increases further towards higher SPS energies. As 
already stated, X^i^t^s) is rather fiat for all energies, so 
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FIG. 5: Freeze-out temperatures {T)i and chemical potentials 
{fJ-B)i of various particle species at -ELab/^ = 30, 80, 158 GeV. 



that S{^b) should also be taken with some uncertainty. 

Fig. 5 shows the freeze-out temperatures and chemical 
potentials for individual particle species at selected ener- 
gies in the CERN-SPS range. The effect of the inhomo- 
geneities is evident. For example, anti-protons are typ- 
ically emitted from regions with lower baryon-chemical 
potential than protons; also, heavy particles are concen- 
trated in "hot spots" while light pions are distributed 
more evenly throughout the decoupling volume etc. The 
temperatures of rare particle species can deviate much 
more from (T) than S{T) from fig. 4 might suggest. This 
is, of course, due to their small weight in eq. (15). Finally, 
we note that the rather high temperatures of the hot 
spots from which the heavy particles emerge might indi- 
cate the need for a better treatment of interactions [13] 
than the simple excluded-volume model employed here. 



SUMMARY AND OUTLOOK 

We have shown that inhomogeneities on the freeze-out 
hypersurface do not average out but reflect in the ^tt ( or 
midrapidity), single-inclusive abundances of various par- 
ticle species. This is due to the non-linear dependence 
of the hadron densities pi{T,fiB) on the local tempera- 
ture and baryon-chemical potential. Consequently, even 
the average probe higher moments of the T and ijb 
distributions. Searching for such inhomogeneities could 
represent a promising observable for a first-order phase 
transition. 

We introduced a simple model for the statistical de- 
scription of chemical freeze-out, extending the commonly 
used grand-canonical ensemble to the case of inhomoge- 
neous decoupling surfaces. If freeze-out occurs shortly 
after a first-order phase transition, hydrodynamic fiow 
and diffusion [4, 34] may not completely wash out inho- 
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mogeneitics produced in the course of the phase transi- 
tion. Given that the number of such "domains" is ex- 
pected to be large, the inclusive distribution of temper- 
ature and baryon-chemical potential on the decoupling 
surface should be approximately Gaussian. 

The model improves the fits of the measured parti- 
cle ratios in the CERN-SPS energy regime -ELabM — 
30 — 160 GeV significantly. Homogeneous freeze-out is 
well outside the 95.4% confidence interval. This suggests 
that at intermediate energies the freeze-out surface is not 
well "stirred". Future studies could perhaps shed more 
light on whether these inhomogeneities can indeed be in- 
terpreted as fingerprints of a first-order phase transition. 
Eventually, one would want to establish more quantita- 
tive relations between the amplitudes of the T, hb inho- 
mogeneities and the properties of the phase transition, 
e.g. its latent heat and interface tension. Furthermore, 
the role of inhomogeneities in the net-strangeness distri- 
bution should be studied. On the other hand, within 
the present model, no statistically significant improve- 
ment over homogeneous freeze-out was observed at lower 
(SPS-20) and higher (RHIC) energies. 

Inhomogeneities could also affect the coalescence prob- 
abilities of (anti-) nucleons to light (anti-) nuclei, which 
are also sensitive to density perturbations [35] . Other sig- 
nals, such as two-particle correlations [8, 36], could also 
be analyzed in this regard. 

To improve the quality of the statistical fits, more data 
on hadron multiplicities would be helpful, in particular 
at the lower end of the CERN-SPS energy spectrum and 
at RHIC. This includes estimates of multiplicities of un- 
stable resonances (p, K*, w, A ...) at chemical freeze- 
out [37]. Data from GSI-FAIR and CERN-LHC will pro- 
vide additional constraints for the evolution of chemical 
freeze-out with energy. 

Note added: While we finished this manuscript, the 

supcrstatistics approach [38] was brought to our atten- 
tion. It considers non-equilibrium systems in stationary 
states with fluctuating intensive quantities, which are de- 
scribed by a superposition of Boltzmann ensembles. Al- 
though our approach emerged from a different physical 
picture, it can nevertheless be viewed as an application of 
"superstatistics" to particle freeze-out in heavy-ion colli- 
sions. 
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